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Abstract We build an ultra-high-speed imaging sys- 

Landt (1967|); Vogel et al (1996); Noack and Vogel|( 1998 

tern based on the background-oriented schlieren (BOS) 

Sankin et al (2008); 

Lauterborn and Vogel (2013 

)). The 

technique in order to capture a laser-induced underwa¬ 

shock waves were applied in various technological fields 

ter shock wave. This BOS technique is able to provide 

such as laser medicine (Hirano et al|(|2002|); Tominaga 

two-dimensional density-gradient field of fluid and re- 

et al (2006); Sankin et al (2005| 2008|)). jTagawa et al 


quires a simple setup. The imaging system consists of 
an ultra-high speed video camera, a laser stroboscope, 
and a patterned background. This system takes images 
every 0.2 /rs. Furthermore, since the density change of 
water disturbed by the shock is exceedingly small, the 
system has high spatial resolution ~ 10 /im/pixel. Us¬ 
ing this BOS system, we examine temporal position of 
a shock wave. The position agrees well with that mea¬ 
sured by conventional shadowgraph, which indicates 
that the high-speed imaging system can successfully 
capture the instantaneous position of the underwater 
shock wave that propagates with the speed of about 
1500 m/s. The local density gradient can be determined 
up to 0(10^ kg/m^), which is confirmed by the gradient 
estimated from the pressure time history measured by 
a hydrophone. 

Keywords Laser-induced underwater shock wave • 
Background-oriented schlieren technique • Ultra 
high-speed recording 


1 Introduction 


(2012) utilized a laser-induced shock wave to generate 


a highly-focused supersonic microjet, which can be ap¬ 


plicable for novel needle-free injection devices (Tagawa 


et al ( 2013[ )). [Turangan (2013) confirmed numerically 
that the laser-induced shock wave plays a key role for 
generating a high-speed microjet. 

The quantification of the shock in experiments is, 
however, challenging due to its ultra-high velocity (about 
1500 m/s) and small change of liquid density. Thus a 
desired measurement system has to have high spatial 
and temporal resolutions. 

In recent years, the background-oriented schlieren 
(BOS) technique has drawn attention as a novel mea- 


surement technique ( 

Dalziel et al (2000); 

Richard and 

Raffel 

(2001); 

Meier 

(2002)). This technique provides 


two-dimensional information of material density with 
quite a simple setup. The technique also has the great 
potential for measuring a wide range of flows. The tech¬ 
nique was applied to a steady supersonic air flow over a 


cone-cylinder model ( 

Venkatakrishnan and Meier ( 

to 

o 

o 

a large-scale imaging outside a laboratory with natu- 

ral backgrounds (Hargather and Settles 

(2010); 

Bar- 


Laser-induced shock waves in water have been exten- generated through a micro-explosion in the air (Suriya- 


sively investigated over several decades (e.g. Bell and 
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narayanan et al (2012); Venkatakrishnan et al (2013)). 
The three-dimensional density field of supersonic air 
flow around an asymmetric body was quantified using 
BOS technique with algebraic reconstruction technique 
(ART) dOta et al| ( pM] )). 

To the best of the authors’ knowledge, an applica¬ 
tion of BOS technique to an underwater shock wave 
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has not been conducted. Here we for the first time uti¬ 
lize BOS technique to measure the underwater shock 
wave. Reminded that the shock has ultra-high velocity 
and is accompanied with little density change in wa¬ 
ter. Thus we construct an ultra-high-speed BOS system 
with an ultra-high-speed video camera, a laser strobo¬ 
scope, and a background. Murphy and Adrian (2011) 
reported that particle image velocimetry (PIV) is appli¬ 
cable to visualize shock waves propagating in transpar¬ 
ent polydimethylsiloxane, which implies that an under¬ 
water shock wave could be measured by BOS technique. 

The structure of this paper is as follows. We describe 
BOS technique in section followed by measurement 
results in section [3l Section |4] concludes our work. 


2 Measurement methods 


2.1 Principle of BOS technique 


Figurel^shows a schematic of BOS technique. The prin¬ 
ciple of the technique is similar to conventional schlieren 
technique which exploits variation of refractive index of 
fluid due to its density gradients (the first-order dif¬ 
ferential value of density; Vp). The variation of re¬ 
fractive index is obtained from two background images 
with random dots. One is a background image with¬ 
out density gradients between the background and the 
image plane (black dash line in Fig. [^. The other is 
a background image with density gradients (red solid 
line in Fig. Q. By using a PIV type cross-correlation 
algorithm, we can obtain background-element displace¬ 
ments. 

Figure shows examples of the obtained images, in 
which thermal plume were captured by a still camera. 
The undisturbed background image and the disturbed 
image are shown in Fig.j^a) and (b), respectively. Fig- 
ure[^c) shows the field of background-element displace¬ 
ments shown as yellow vectors. 

These background-element displacements are related 


to density gradients as (Venkatakrishnan and Meier 

( 2 ^); 


V = 


L 


Zd+AZo 

—dz 


Zb ^0 J Zd-^Zd 


n = Kp -f-1, 


( 1 ) 

( 2 ) 


where v is the background-element displacement in the 
vertical axis {y direction, see Fig.[^, Zb is the distance 
from the background to the schlieren object, / is the 
focal length, is the distance from the background to 
the lens, AZd is one half of the thickness for schlieren 
object, n is the refractive index, z is the line-of-sight 
direction, p is the density, and K is the Gladstone-Dale 


constant. Note that the background-element displace¬ 
ments obtained using BOS technique (Eq.[^ are related 
to the integrated value of the density gradients. 


Background 



Fig. 1 Schematic of BOS technique. A black dash line shows 
a ray from a background image without density gradients 
between the background and the image plane. A red solid 
line shows a ray with density gradients 



Fig. 2 Background images (a) without density gradients and 
(b) with density gradients due to a flame, (c) The image shows 
the difference between (a) and (b) using a PIV-type cross¬ 
correlation algorithm. Yellow vectors show the local displace¬ 
ments of the background 


2.2 Experimental setup and conditions 

Eigurej^ shows our experimental setup. Our BOS setup 
consists of a background with random dot pattern, a 
schlieren object (a density gradient), an ultra-high-speed 
camera, and a back light. A 532 nm, 6 ns, laser pulse 
(Nd: YAGlaser Nano S PIV, Litron Lasers ltd., UK) 
is focused through a microscope objective lens (SLM- 
PLN20x, magnihcation: 20x, N.A.: 0.25, Olympus co., 
Japan) to a point inside a water-hhed pool (450 x 300 x 
300 mm), inside which a shock wave emerges. We record 
images using an ultra-high-speed video camera (Kirana, 
Specialized Imaging co., UK) with up to 5 x 10^ frames 
per second and 924x768 pixel array. The experiment 
requires exceedingly short exposure time since the flow 
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is non-stationary. Hence we utilize a laser stroboscope 
(CAVILUX Smart, CAVITAR co., Finland) with a pulse 
width of 20 ns as an illumination source. Its repeti¬ 
tion rate is also up to 5 MHz. A digital delay generator 
(Model 575, BNC co., USA) synchronizes the camera, 
the laser, and the stroboscope. 

We obtain background-element displacements using 
a PIV-type cross-correlation algorithm. We use 50% 
overlap with the Fast Fourier Transform (FFT) multi 
pass interrogation (from 32 pixel to 4 pixel). 

For validating the high-speed BOS system, we con¬ 
duct two kinds of experiments. 

First, we compare temporal positions of a shock 
wave measured by BOS technique with that by shadow¬ 
graph, a common conventional method to visualize an 


underwater shock wave (Bell and Landt (1967); Sankin 
et al (2008); Brujan et al| ( |2008| )). Shadowgraph tech¬ 


nique consists of a simple optical setup; a schlieren ob¬ 
ject (a density gradient), a camera, and a back light. 
When we take shadowgraph images with our setup, we 
just remove the background from the high-speed BOS 
system. Note that shadowgraph records local bright¬ 
ness that corresponds to second derivative of the fluid 
density (Settles (2001)). Since the displacement 
obtained from BOS technique is related to the density 
gradient Vp, we differentiate the density gradient for 
obtaining second derivative of the fluid density V^p. In 
order to extract the mean shock position we approx¬ 
imate the shape of the shock as a sphere (a circle in 
two-dimensional projection). 

Second, we verify the displacements obtained from 
our BOS technique by comparing them with those es¬ 
timated from hydrophones. We use two hydrophones 
(Muller-Platte Needle Probe, Muller ins., Germany) and 
the high-speed BOS system simultaneously. Each hy¬ 
drophone is arranged at right angles to the direction of 


Fig. 4 Time history of shock pressures measured by two hy¬ 
drophones arranged at right angles to the direction of a laser 
beam at a stand-off distance of 5 mm. Input laser energy is 
1.9 mJ 


Shock pressure 





◄- 

Displacement related 
to pressure gradient 


(b) Reconstructed 
shock front 



Fig. 5 Procedure of converting shock pressure to displace¬ 
ments. (a) Conversion of shock pressure measured by a hy¬ 
drophone to displacements related to pressure (density) gra¬ 
dients. (b) Reconstructed displacement field with assuming 
an axisymmetric shock structure, (c) Projection of the recon¬ 
structed displacement field toward 2 :-direction 


a laser beam at a stand-off distance of about 5 mm 
and is connected to an oscilloscope (ViewGo H DS- 
5554A, Iwatsu co., Japan). Figureshows typical hy¬ 
drophone measurements. Both shock pressures are al¬ 
most the same, which indicates that the shock wave 
has an axisymmetric structure. The small difference is 
due to different positions of hydrophones; a stand-off 
distance of the hydrophone #1 is 5.14 mm and the hy¬ 
drophone ^2 is 5.27 mm. 

A procedure for estimating displacements from hy¬ 
drophone measurements is as follows: We first calcu¬ 
late pressure gradient from pressure time history (Fig. 
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I^a)). We obtain the refractive index using Gladstone- 
Dale equation (Eq. and Tait equation; 

P±^=(pX ( 3 ) 

Po + B (poj ’ ^ ^ 

where po is the hydrostatic pressure, B and m are con¬ 
stants based on physical properties. In the case of wa¬ 
ter, B is 314 MPa and m is 7. The calculated varia¬ 
tion of refractive index is then related to displacements 
via Eq. We then reconstruct the two-dimensional dis¬ 
placement field on x-z plane from hydrophone results as 
shown in Eig. [^b) since a hydrophone provides point- 
wise pressure information while the displacement ob¬ 
tained from BOS technique is a value integrated along 
z-axis as shown in Eig. i c) (cf. Fig. 0 . Reminded that 
we here assume an axisymmetric structure of a laser- 
induced shock wave, which is verified from hydrophone 
measurements in Eig. Einally, we obtain the esti¬ 
mated displacements by integrating the reconstructed 
displacement fields along z-axis. 

We set two levels of laser energy: 1.9 mJ and 3.0 
mJ. In addition, we change spatial resolutions of im¬ 
ages, which depend on the focal length / and the dis¬ 
tance from the background to the lens Zb- With higher 
spatial resolution of 8.4 /im/pixel, / and Zb are respec¬ 
tively 360 mm and 120 mm, while with lower spatial 
resolution of 14.5 /im/pixel, / and Zb are respectively 
250 mm and 120 mm. 


3 Results and Discussion 


3.1 Position of an underwater shock wave visualized 
by BOS technique and shadowgraph 


We show spatiotemporal evolution of a shock wave mea¬ 
sured by BOS technique and shadowgraph in Eig. 
Original images obtained from the high-speed BOS sys¬ 
tem are displayed in Eigurej^a). Eigurej^b) shows the 
vector field of displacements of background-elements, 
which correspond to the density gradients. The regions 
of non-zero density gradients (non-zero vectors) radi¬ 
ally propagates from a laser-induced bubble. Eor clear 
presentation, the magnitudes of the displacements are 
shown in Eig.j^c). Red color shows the largest displace¬ 
ment. We observe spherical-symmetric distribution of 
the shock as reported in many previous researches (e.g. 


Vogel et al (1996); [Petkovsek et al (2005); Brujan et al 
(2008)). Eigure |6 M) shows the second-order differential 
value of density (V^p) obtained by differentiating the 
density gradient filed (Vp). Eor comparison, a shock 
wave in the same condition visualized by shadowgraph 
is shown in Eig. [^e). A dark streak pattern due to a 
shock wave shows the second-order differential value of 


density (V^p). Both Eig. |^d) and (e) look almost the 
same, showing sphere-like shapes of the shocks. 

We here quantitatively compare shock positions by 
both methods. We plot both shock radii, namely mean 
shock positions, as a function of time in Eig. The 
temporal evolution of the shock position from the BOS 
system agrees with that by shadowgraph, within mea¬ 
surement error ±120 /im. 

Eor a spherical shock wave, the temporal evolution 
can be described as ( |Dewey| ( 2QQ1| )), 

Rs — Cl ± C 20 jot ± Cs ln(l ± dot) ± C74y/ln(l ± ngt), (4) 


where Rs is the htted shock position, Ci, C 3 , and C 4 
are the htting constants, C 2 is set to 1 , ao is the speed 
of sound in water (ao=1483 m/s), and t is an elapsed 
time after the laser being hred. It nicely hts the data, 
indicating that the shock visualization is valid. Thus 
the high-speed BOS system is able to detect the in¬ 
stantaneous position of the non-zero density gradient, 
i.e. an underwater shock wave. 


3.2 Comparison of BOS technique with hydrophone 
measurement 

In this section we compare the displacements obtained 
from BOS technique with those calculated from hy¬ 
drophone measurements. The displacement hleds of x- 
components (u) for two different laser energies are shown 
in the upper side of Eigs. |^and|^ The spatial resolu¬ 
tion is 8.4 /im/pixel. The shape of the shocks are almost 
spherical for ah the laser energies. Eor quantitative dis¬ 
cussion, the displacements u are shown in the lower 
side of Eigs. and where the displacements obtained 
from BOS technique and a hydrophone are respectively 
presented as a black line and a red line. Shock pres¬ 
sure measured by a hydrophone is also displayed as a 
blue-chain line for reference. 

Eor lower laser energy (1.9 mJ), the maximum dis¬ 
placement obtained from BOS technique fairly agree 
with that from hydrophone measurements. In this par¬ 
ticular case the maximum density gradient is 6 x 10 ^ 
kg/m^ based on peak pressure of 1 MPa and shock 
width of 60 /im. 

As the laser energy increases, however, the maxi¬ 
mum displacement obtained by BOS technique is con¬ 
siderably smaller than the displacement calculated from 
hydrophone measurements. The ratios of the maximum 
displacement from BOS technique to that from hydrophones 
are 87% and 26% for laser energy 1.9 mJ and 3.0 mJ, 
respectively. It indicates that in the cases for higher 
laser energy, the largest density gradient exceeds a mea¬ 
surable limit of the high-speed BOS system. Eor input 
laser energy 3.0 mJ, the maximum density gradients are 
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Fig. 6 The spatiotemporal evolution of the laser-induced shock wave after the laser pulse beam being fired. The laser pulse 
beam is illuminated from the top of these images, (a)-(e) are shown at the same timing. We mask the region of bubble with 
brown elliptic area in (b)-(d). (a) Original images obtained from the ultra-high-speed recording system, (b) Displacement field 
obtained by BOS technique. The displacements are shown as yellow vectors. The displacement corresponds to the density 
gradient (Vp). (c) Magnitude of the density gradient field (|Vp|). (d) The second-order differential value of density (V^p) 
calculated from the displacement values (b). (e) The shock wave visualized by shadowgraph. A dark streak pattern due to 
shock wave shows the second-order differential value of density (V^p) 
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Time, jjs 


Fig. 7 Propagation of the shock front plotted as a function 
of time . The s olid line presents the fitting curve proposed by 
|Dewey| ( |2001| ) 


1.4x10^ kg/m^. Therefore current high-speed BOS sys¬ 
tem can be applicable to a flow field, density gradient 
of which is of the order of 10^ kg/m^. 

In order to verify the effect of spatial resolution 
of our BOS system, we measure the underwater shock 
wave with lower resolution (14.5 /im/pixel) in the con¬ 
dition of 1.9 mJ laser energy. Measured displacements 
are shown in Fig. The maximum displacement for 
14.5 /im/pixel is 46% of that from a hydrophone, which 
is obviously less than the value obtained with higher 
spatial resolution (8.4 /im/pixel). We thus confirm that 
the applicable range of the BOS system depends on its 
spatial resolution as pointed out by |Goldhahn and Se- 
(120071). 


ume 


4 Conclusion 

We have applied background-oriented schlieren (BOS) 
technique to an underwater shock wave. Challenges in 
the measurements are ultra-high velocity of the shock 
wave in water and small density fluctuation of water. 
We have constructed an ultra-high speed BOS system 
that is able to take images every 0.2 /is with high-spatial 
resolution ^ 10 /im/pixel. For its validation, we have 
measured the shock by shadowgraph or hydrophones 
simultaneously with the BOS system. 

We have found that the spatiotemporal position of 
the shock detected by BOS technique has agreed with 
that by shadowgraph. It indicates that the high-speed 
BOS system is able to detect instantaneous position 
of the shock wave that changes local density of water. 
In addition, in order to estimate a measurable range 
for density gradient we have compared the displace¬ 



Distance, mm 


Fig. 8 Displacements obtained from BOS technique and 
those calculated from hydrophone measurements with 1.9 mJ 
of laser energy. Upper figure presents displacement fields of 
the x-direction (u) obtained from BOS technique with spatial 
resolution of 8.4 jum/pixel. Lower figure compares displace¬ 
ments. The black line represents displacements obtained from 
BOS technique; the red line shows displacements calculated 
from hydrophone measurements; the blue dashed line is shock 
pressure measured by a hydrophone 


ments obtained from BOS technique with those from 
hydrophones. We have varied input laser energy and 
spatial resolution of the BOS system. The measurable 
range increases with an increase of spatial resolution of 
the BOS system. It has been found that with current 
high-speed BOS system the density gradient of the or¬ 
der of 10^ kg/m^ can be reasonably measured. Spatial 
resolution higher than 8.4 /im/pixel might enable us to 
measure a flow field with even higher density gradients. 
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Distance, mm 


Fig. 9 Upper panel: Magnitude of the density gradient field 
presented in the same manner for Fig. [^c). Lower panel: 
Displacements obtained from BOS technique and those cal¬ 
culated from hydrophone measurements. The laser energy is 
3.0 mJ. Lines in this figure represent the same quantities as 
those in Fig. 
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